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LONG-TERM  GOALS 

The  long-term  goals  of  this  research  project  are  to  advance  data  assimilation  (DA)  methods  for  coastal 
ocean  circulation  models.  After  several  years  of  research  focused  on  various  aspects  of  advanced, 
variational  DA  with  nonlinear  coastal  ocean  models,  we  have  come  close  to  development  of  an 
optimal,  versatile,  and  relocatable  DA  system  based  on  a  primitive  equation  model  with  a  turbulence 
submodel.  The  system  can  be  used  efficiently  both  for  operational  needs  (forecasting,  search  and 
rescue,  environmental  response)  and  for  fundamental  studies  of  coastal  ocean  dynamics.  The  modeling 
system  will  have  the  capability  to  assimilate  time-series  measurements  from  moorings  and  coastally- 
based  high  frequency  (HF)  standard-  and  long-range  radars,  satellite  information  (SSH,  SST),  and 
hydrographic  survey  data  (including  observations  from  underwater  gliders). 

OBJECTIVES 

Scientific  and  technological  objectives  of  this  project  include: 

Development  of  a  data  assimilation  system  that  utilizes  the  advanced  variational  representer 
method  with  a  nonlinear  three-dimensional  model  of  ocean  circulation  over  shelf  and  in  the 
adjacent  interior  ocean  energized  by  coastal  flows  (a  coastal  transition  zone,  CTZ). 

Testing  this  data  assimilation  system  with  observations  available  along  the  Oregon  coast. 
Understanding  the  impact  of  assimilating  satellite  observations  (alongtrack  SSH,  SST  composites) 
and  in-situ  observations  (HF  radar  surface  velocities,  moored  velocities,  glider  hydrographic 
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transects)  to  improve  accuracy  of  prediction  of  the  oceanic  currents  and  hydrographic  conditions, 
including  short  term  (3-6  day)  forecasts. 

Understanding  co-variability  of  shelf  and  interior  ocean  processes,  using  assimilation  results  and 
computations  of  representers  (elements  of  the  time-  and  space-dependent  model  error  covariance). 

APPROACH 

The  proposed  research  involves  a  systematic  continuation  of  work  in  progress  that  has  included 
assimilation  of  current  measurements  from  moored  instruments,  surface  velocities  from  an  array  of  HF 
radars  deployed  along  the  Oregon  coast,  and  satellite  SSH  and  SST.  Additional  data  types,  including 
in-situ  sections  of  hydrographic  and  turbulence  observations,  have  been  used  to  verify  results  of  data 
assimilation.  Studies  of  the  data  assimilation  methodologies  have  been  focused  on  the  effectiveness  of 
various  methods,  from  simple  (albeit  suboptimal)  “optimal  interpolation”  to  more  advanced  variational 
methods.  These  methods  have  been  utilized  with  models  of  increased  complexity  (linear  models,  a 
nonlinear  shallow  water  model,  and,  presently,  a  three-dimensional  nonlinear  model  of  stratified  flows 
with  a  turbulence  submodel).  The  data  assimilation  results  have  been  reported  by  Scott  et  al.,  (2000), 
Oke  et  al.  (2002a),  and  Kurapov  et  al.  (1999,  2002,  2003,  2005a,  2005b,  2005c,  2007,  2009,  2010b). 
These  data  assimilation  studies  have  been  complemented  by  the  dynamical  analyses  of  wind-driven 
upwelling  (e.g.,  Oke  et  al.,  2002b,  2002c,  Kurapov  et  al.,  2005c,  Springer  et  al.,  2009),  jets  in  the  CTZ 
(Koch  et  al.,  2010),  and  internal  tides  in  combination  coastal  upwelling  (Kurapov  et  al.,  2003,  2010a). 

Our  present  approach  is  based  on  the  use  of  the  variational  representer  method  (Chua  and  Bennett, 
2001,  Bennett,  2002)  with  the  nonlinear  Regional  Ocean  Modeling  System  (ROMS),  a  free-surface, 
hydrostatic,  primitive-equation  model  featuring  terrain-following  coordinates  and  advanced  numerics 
(Shchepetkin  and  McWilliams,  2005).  The  variational  method  with  a  nonlinear  model  requires 
implementation  of  a  companion  tangent  linear  model  and  its  adjoint  counterpart.  We  have  developed 
and  tested  our  own  version  of  the  tangent  linear  and  adjoint  codes,  AVRORA  (Advanced  Variational 
Regional  Ocean  Representer  Analyzer).  The  dynamics  of  these  stand-alone  codes  are  consistent  with 
the  nonlinear  ROMS.  Availability  of  our  own  tangent  linear  and  adjoint  codes  has  allowed  us 
considerable  flexibility  in  the  choice  of  the  critical  details  of  our  data  assimilation  system,  including  (i) 
ways  in  which  observations  are  interpreted  by  the  model  (data  functionals)  and  (ii)  nontrivial 
hypotheses  about  error  covariances  for  the  corrected  initial  conditions  and  atmospheric  forcing. 
Although  tangent  linear  &  adjoint  ROMS  codes  have  been  distributed  to  the  research  community, 
those  codes  do  not  allow  us  similar  flexibility.  In  particular,  some  of  the  data  functionals  that  we  have 
tested  (e.g.,  satellite  alongtrack  SSH  slope,  HF  radar  surface  current  radial  component  data,  any  time- 
averaged  data)  are  not  built  in  ROMS. 

In  variational  data  assimilation,  the  improved  ocean  state  estimate  is  found  by  solving  a  least-squares 
problem.  The  optimal  ocean  state  minimizes  a  cost  function  defined  as  a  sum  of  quadratic  penalty 
tenns  on  errors  in  model  inputs  (generally  including  initial  conditions,  boundary  conditions,  forcing, 
also  possibly  time-  and  space-distributed  errors  in  model  equations)  and  data-model  misfits,  all 
integrated  over  space  and  a  specified  time  interval.  Fonnulation  of  each  penalty  tenn  requires 
specification  of  a  corresponding  error  covariance.  In  particular,  the  error  covariances  for  the  model 
inputs  can  provide  smoothing  (filtering)  of  the  model  correction  and  also  possibly  dynamical  balances 
in  the  corrected  fields.  To  find  the  optimum  solution,  we  utilize  the  indirect  representer  algorithm 
(Chua  and  Bennett,  2002)  with  preconditioning  (Egbert  et  al.,  1994)  such  that  the  method  is  applicable 
with  large  data  sets.  As  for  any  variational  method  applied  to  a  nonlinear  dynamical  model,  the 
representer  method  requires  repeated  solution  of  the  corresponding  tangent  linear  and  adjoint  systems. 
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Collaborators  in  this  project  include  Dr.  Kurapov  (lead  PI),  the  developer  of  the  AVRORA-ROMS 
assimilation  system  and  the  lead  PI.  The  co-PIs  Drs.  J.  S.  Allen  and  G.  Egbert  have  provide  valuable 
expertise  and  advice  on  coastal  ocean  dynamics  and  variational  data  assimilation.  Dr.  Peng  Yu  joined 
this  project  in  February  2009  to  help  accomplish  massive  computations  and  analyses. 

WORK  COMPLETED 

Our  studies  have  been  focused  on  an  area  centered  on  the  Oregon  shelf  and  CTZ  extending  between 
41-4 7N  and  124-129W  (Figure  1).  In  preparation  to  assimilation  tests,  a  3  km  resolution  ROMS  model 
has  been  set.  Boundary  information  has  been  provided  from  the  9-km  resolution  Navy  Coastal  Ocean 
Model  of  the  California  Current  System  (NCOM-CCS)  run  at  NRL  (Shulman  et  ah,  2009).  Model 
studies  without  data  assimilation  were  leveraged  by  efforts  on  the  recent  NOPP  and  GLOBEC 
modeling  projects  (Springer  et  ah,  2009;  Koch  et  ah,  2009).  Assimilation  tests  have  been  perfonned  at 
the  6-km  horizontal  resolution,  with  focus  on  the  impact  of  alongtrack  SSH  altimetry,  HF  radar  surface 
currents,  and  satellite  SST  on  surface  and  subsurface  oceanic  fields. 


Ion,  W 

Figure  1.  Oregon  model  domain  with  locations  of  observations  available  recently  (dots  -  altimetry 
tracks,  gray  area  -  HF  radar  surface  currents,  blue  lines  -  glider,  star  -  mooring). 

Assimilation  tests  of  the  AVRORA-ROMS  system  with  actual  observations  have  been  performed  in  a 
series  of  relatively  short  windows  (1,  3,  or  6  days).  In  each  window,  the  variational  AVRORA  system 
is  run  to  find  correction  to  the  initial  conditions  at  the  beginning  of  the  window.  Then,  the  nonlinear 
ROMS  is  run  for  the  period  twice  as  long  as  the  assimilation  interval,  to  provide  the  analysis  solution 
(in  the  assimilation  window)  and  the  forecast  (to  be  used  as  the  prior  model  in  the  next  assimilation 
window). 

The  study  focused  on  assimilation  of  alongtrack  altimetry  and  multi-satellite  blended  SST  (SST  data 
provided  by  D.  Foley  from  NOAA)  has  been  completed  for  the  study  period  of  2005,  when  data  from 
three  satellites  (Topex,  Jason,  Envisat)  were  available  (see  Figure  1).  The  manuscript  describing  this 
study  was  submitted  to  J.  Geophys.  Res.  (Kurapov  et  al.,  2010b)  and  is  currently  at  the  revision  stage. 
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Dr.  P.  Yu’s  studies  have  been  focused  on  assimilation  of  HF  radar  surface  currents  (data  provided  by 
P.  M.  Kosro,  OSU).  Tests  have  been  completed  for  summer  2008  and  the  manuscript  is  in  preparation. 

The  work  assimilating  actual  observations  has  been  complemented  by  the  analyses  of  the  representers 
as  functions  that  show  zones  of  influence  of  assimilated  observations.  In  particular,  in  these  analyses 
we  have  looked  at  the  effect  of  the  initial  condition  error  covariance  (balanced  vs.  unbalanced)  and  also 
at  the  propagation  of  the  coastally  trapped  waves  that  would  carry  assimilated  infonnation  from  south 
to  north  in  a  realistic  coastal  environment. 

Data  assimilation  studies  have  been  complemented  by  new  analyses  of  the  nonlinear  ocean  dynamics 
off  Oregon.  Using  leveraging  by  this  and  another  projects,  we  have  completed  the  study  of  the  near¬ 
surface  dynamics  of  the  energetic  jet  separated  in  the  coastal  transition  zone  off  Cape  Blanco  (south 
Oregon)  (Koch  et  ah,  2010).  Kurapov  et  al.  (2010a)  have  completed  the  study  on  the  intermittency  of 
the  internal  tide  and  its  influences  on  the  wind-driven  currents  on  the  shelf. 

RESULTS 

Tests  of  the  AVRORA-ROMS  assimilation  system  with  a  6-km  resolution  coastal  ocean  model  and 
actual  observations  of  HF  radar  surface  velocities,  satellite  SSH  and  SST  have  demonstrated  feasibility 
of  the  application  of  the  variational  representer-based  method  in  a  real-time  coastal  ocean  forecast 
system.  After  a  preconditioner  is  built  (requiring  computation  of  100-200  representers,  running  adjoint 
and  tangent  linear  jobs  in  parallel),  the  convergence  rates  are  good,  requiring  only  5-15  iterations  of  the 
tangent  linear  and  adjoint  AVRORA  models  in  each  assimilation  window  to  obtain  sensible 
improvement  of  the  initial  conditions.  With  further  improvements  in  parallelization  and  minimization 
algorithm,  assimilation  can  be  approached  at  better  (2-3  km)  resolution. 

In  the  study  assimilating  satellite  alongtrack  SSH  (Kurapov  et  al.,  2010b),  we  choose  to  fit  the  model 
SSH  slope  to  the  observed  slope.  This  allows  avoiding  a  problem  of  matching  mean  levels  in  the  model 
and  observations.  As  the  effect  of  SSH  slope  assimilation,  surface  geostrophic  currents  are  improved. 

In  addition,  as  verified  by  comparison  with  the  unassimilated  SST,  the  geometry  of  the  upwelling  SST 
front  is  improved  dramatically  (Figure  2).  The  SSH  assimilation  positively  affects  subsurface  fields  as 
verified  against  (limited)  hydrographic  section  data.  Later  in  summer  2005,  the  outer  front  is 
developed,  extending  parallel  to  the  coast  at  a  distance  of  approximately  200  km.  This  front  is  not 
reproduced  in  the  limited  resolution  free-run  model  (i.e.,  the  model  that  does  not  assimilate  data)  and  is 
not  always  apparent  in  the  SST  maps,  since  it  is  rather  associated  with  strong  subsurface  horizontal 
density  gradient. 
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Figure  2.  Assimilation  of  satellite  alongtrack  SSH  data  in  the  coastal  transition  zone  off  Oregon 
qualitatively  improves  the  geometry  of  the  SST  front:  (left)  prior,  free  run  6-km  resolution  model, 
(center)  model  constrained  by  assimilation  of  SSH  altimetry  beginning  June  2005,  (right) 
verification  GOES  SST  daily  composite  (all  images  are  for  25  September  2005). 
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Figure  3.  Assimilation  of  HF  radar  surface  currents  helps  to  improve  prediction  of  not  only 
currents,  but  also  SST:  (a)  GOES  SST  (degrees  C)  and  HF  radar  mapped  surface  velocities,  day 
205,  2008,  (b)  prior  model  SST  and  currents  (shown  at  the  same  locations  as  observed  currents), 
and  (c)  forecast  (SST  and  surface  currents),  after  initial  conditions  on  day  201  were  corrected 
assimilating  HF  radar  currents  for  days  201-203  [P.  Yu]. 
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Common  for  most  assimilation  methods,  the  series  of  analyses  solutions  exhibits  instantaneous  finite 
amplitude  steps  every  time  the  data  assimilation  correction  to  initial  conditions  is  made.  This 
correction  can  be  interpreted  as  the  impulse  force  added  on  the  r.h.s.  of  equations.  In  other  words,  the 
DA  solutions  are  generally  not  dynamically  balanced  in  the  time-averaged  sense.  To  evaluate  the 
relative  magnitude  of  the  correction  in  heat  content  resulting  from  assimilation  of  the  alongtrack  SSH 
slope  or  the  slope  and  SST  in  combination  we  have  perfonned  the  analysis  of  the  volume  integrated 
heat  budget  near  coast  (including  shelf  and  a  part  of  the  interior  ocean).  We  find  that  the  DA  correction 
tenn  is  approximately  of  the  same  magnitude  as  the  net  horizontal  advection  of  heat  and  atmospheric 
heat  flux.  Assimilation  of  SST  in  addition  to  SSH  not  only  improves  the  SST  values,  but  also  helps  to 
reduce  the  negative  bias  in  the  DA  correction  term. 

Assimilation  of  HF  radar  surface  velocities  (daily  averaged  mapped  velocity  fields)  in  a  series  of  3-day 
time  windows  helps  to  improve  accuracy  of  both  analysis  and  3-day  forecast  surface  velocity  fields. 
Velocity  assimilation  also  helps  to  improve  the  SST  structure  (Figure  3),  as  confirmed  by  statistical 
comparisons  with  the  unassimilated  daily  SST  fields.  Each  HF  radar  provides  information  on  a 
component  of  the  surface  velocity  in  the  direction  of  the  radar.  We  have  run  tests  assimilating  the 
radial  component  data  and  results  were  similar  to  the  case  assimilating  mapped  (u,v)  data.  This 
experiment  showed  that  the  variational  assimilation  system  can  be  used  as  a  new  tool  for  interpolation 
and  mapping  HF  radar  radial  component  data  from  a  number  of  radars,  to  obtain  maps  of  orthogonal 
velocity  components  at  each  location.  The  advantage  of  such  a  mapping  tool  is  that  it  involves 
dynamically  based  scales  in  the  interpolation  functions,  provides  opportunity  to  synthesize  the  HF 
radar  data  with  other  data  types,  and  yields  accurate  forecasts  of  surface  currents  and  SST. 

Success  in  assimilation  of  multi-satellite  blended  SST  products  can  be  limited  because  the  noise  in  the 
data  from  each  satellite  often  contaminates  the  blended  product.  This  is  in  particularly  true  of  GOES 
geostationary  satellite  SST.  To  alleviate  this  problem,  we  have  run  tests  assimilating  original  hourly 
GOES  data  in  the  variational  system.  Time-  and  space-interpolation  of  these  sparse  and  noisy  fields, 
with  projection  of  the  correction  to  the  initial  conditions  at  the  beginning  of  the  3-day  time  window, 
allows  removing  high-frequency  noise  in  the  hourly  GOES  data  and  filling  gaps. 

Analyses  of  the  results  of  the  adjoint  solution  confirmed  that  the  zones  of  influence  of  the  data  are 
affected  by  the  dynamical  processes,  including  advection  and  coastally  trapped  waves  (Figure  4). 
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Figure  4.  The  SSH  component  of  the  adjoint  solution  for  the  observation  of  the  alongtrack  SSH 
slope  at  the  location  shown  as  square,  showing  that  the  error  in  the  observed  model  variable  at  later 
time  (t=3  d)  co-varies  with  the  initial  error  in  the  SSH  south  of  the  observational  location,  as  an 

effect  of  coastally  trapped  waves. 


IMPACT/APPLICATIONS 

Implementation  of  the  variational  data  assimilation  system  to  the  Oregon  shelf  and  adjacent  interior 
ocean  has  demonstrated  value  of  satellite  alongtrack  SSH,  SST,  and  HF  radar  observations  for  accurate 
multivariate  prediction  and  forecasting  in  the  coastal  ocean.  Understanding  evolution  of  the  model 
error,  detected  by  the  tangent  linear  model  and  its  adjoint,  helps  explain  covariability  of  different 
processes  in  the  ocean  (surface  and  subsurface,  coastal  and  interior  ocean). 

TRANSITIONS 

The  data  assimilation  system  developed  in  this  project  has  become  a  part  of  the  pilot  real-time  Oregon 
coastal  ocean  prediction  system  that  provides  daily  updates  of  3-day  forecasts  in  the  area  centered  on 
the  Oregon  coast  (the  real-time  work  has  been  leveraged  by  the  NOAA  grants).  Since  August  2009,  the 
data  assimilation  system  works  in  the  real-time  operational  mode,  assimilating  HF  radar  surface 
current  and  GOES  hourly  SST.  This  prediction  system  is  a  part  of  the  regional  integrated  ocean 
observing  system  (IOOS)  supported  by  NOAA  via  the  NANOOS  regional  association. 
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RELATED  PROJECTS 


ONR,  “Combining  sequential  and  variational  data  assimilation”.  In  this  new  project,  we  collaborate 
with  Dr.  Peter  Oke  (SCIRO,  Australia)  to  provide  comparative  tests  of  optimal  interpolation  (01)  and 
variational  data  assimilation  approaches.  The  study  will  involve  transition  of  AVRORA-ROMS  to  our 
Australian  partners. 

NSF,  “Modeling,  assimilation,  and  analysis  of  the  shelf-  open  ocean  exchange  off  Oregon.”  The  data 
assimilation  system  utilized  and  advanced  in  this  project  will  be  used  to  study  effect  of  surface 
assimilation  on  subsurface  fields  and  ultimately  effects  of  climate  variability  on  coastal  ocean 
circulation. 

NOAA-IOOS,  “Enhancing  the  Pacific  Northwest  Regional  Coastal  Observing  System  (RCOOS)  of 
Northwestern  American  Network  of  Ocean  Observing  Systems  (NANOOS)”.  Data  assimilation  system 
improvements  will  be  incorporated  in  the  pilot  real-time  ocean  forecast  system  run  operationally  off 
Oregon  with  assimilation  of  satellite  alongtrack  SSH,  GOES  SST,  and  HF  radar  surface  currents. 
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